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New ab initio potential energy surfaces for the 2 II ground electronic state of the Ar-SH complex are 
presented, calculated at the RCCSD(T)/aug-cc-pV5Z level. Weakly bound rotation-vibration levels 
are calculated using coupled-channel methods that properly account for the coupling between the 
two electronic states. The resulting wavefunctions are analysed and a new adiabatic approximation 
including spin-orbit coupling is proposed. The ground-state wavefunctions are combined with those 
obtained for the excited 2 E + state [Phys. Chem. Chem. Phys. 6, 5463 (2004)] to produce transition 
dipole moments. Modelling the transition intensities as a combination of these dipole moments and 
calculated lifetime values [J. Chem. Phys. 109, 170 (1998)] leads to a good representation of the 
experimental fluorescence excitation spectrum [J. Chem. Phys. 98, 4301 (1993)]. 



I. INTRODUCTION 

Van der Waals complexes containing open-shell species 
are of great current interest. In particular, complexes 
containing atoms or molecules with orbital angular mo- 
mentum necessarily involve multiple electronic statesi^ 
They provide a test-bed for studying electronically non- 
adiabatic effects, which are important in the theory 
of reaction dynamics In addition, the observation 
of pre-reactive Van der Waals complexes trapped in 
bound levels§i&2iiii can shed light on intermolecular 
forces in the entrance and exit channels of chemical 
reactionsiii*i2iiii4 The form of these shallow, long- 
range wells can be important in determining reaction 
outcome a 15 i 16 and transition-state geometriesi 17 i 18 

In this paper, we consider the complex consisting of 
an open-shell SH radical and an Ar atom. New ab ini- 
tio potential energy surfaces (PESs) for the X 2 n state 
are presented, and used in calculations of the bound 
rotation-vibration levels. We discuss the possibility of 
employing an approximate approach in the bound-state 
calculations, using a single adiabatic PES rather than 
the two surfaces used in the standard methodi Finally, 
the bound-state energies and wavefunctions are used to 
simulate the vibrationally resolved electronic spectrum. 

The Ar-SH cluster was first detected experimentally 
by Miller and co-workers^ using laser-induced fluores- 
cence excitation spectroscopy. Subsequently this group 
developed empirical potential energy surfaces for the 
complex, for both the A stated and the X stated by 
fitting model functions to reproduce laser-induced fluo- 
rescence (LIF) results. The region of the A state PES 
corresponding to the Ar-S-H configuration (Jacobi an- 
gles between ~ 90 and 180 degrees) was determined only 



approximately, because the fluorescence experiments did 
not probe this zone. More recently, Hirst et al? 2 have 
presented a PES for the A state based on ab-initio cal- 
culations at the RCCSD(T) level with the aug-cc-pV5Z 
basis set. This surface was used to predict bound vibra- 
tional levels in the Ar-S-H configuration 22 which have 
not, to our knowledge, been observed in experiment so 
far. A possible reason why these levels have eluded de- 
tection is discussed in section IV of this paper. 

Sumiyoshi et alm^ have recorded high-resolution spec- 
tra for Ar-SH in the ground electronic state us- 
ing Fourier-transform microwave (FTMW) spectroscopy. 
These authors also produced PESs for the X state based 
on fitting a function to reproduce their experimental 
results^ and these surfaces were later improved with 
the aid of some ab initio results^ Most recently, re- 
sults from microwave-millimctcr-wave double-resonance 
spectroscopy 2 ^ were employed to determine new three- 
dimensional PESs for the X stated 

The family of weakly bound clusters containing a rare 
gas atom and either the OH or SH radical have been 
reviewed by Carter et al£L in 2000 and by Heaven 2 ^ in 
2005. 

The structure of the present paper is as follows. In 
Section II we present new PESs for Ar-SH ( 2 n) based 
entirely on ab initio calculations at the RCCSD(T) level 
with an aug-cc-pV5Z basis set. In Section III we de- 
scribe bound rotation-vibration level calculations using 
these surfaces. We also investigate the wavefunctions and 
introduce a new adiabatic approximation for the bound 
states, including spin-orbit coupling. In Section IV the 
results are combined with those of a previous study of 
the A 2 E + state, in order to produce a high-quality sim- 
ulation of the vibrationally resolved electronic spectrum. 



II. POTENTIAL ENERGY SURFACES 

The geometry of the complex is specified in terms of 
body-fixed Jacobi coordinates r, R and 9. R is the length 
of the vector R which links the center of mass of the 
SH fragment to the Ar nucleus. The vector r links the 
S nucleus to the H nucleus: its modulus r is the SH 
bond length. The angle between R and r is 9, so that 
9 = corresponds to a linear Ar-H-S configuration. For 
this work the bond length r was held constant at the 
experimentally determined equilibrium value of 1.3409 
Aj22 which is justified because the vibrational motion of 
the diatom is very weakly coupled to the relatively low- 
frequency Van der Waals modes of interest. 

Energies were calculated for a regular grid of ge- 
ometry points using the MOLPRO quantum chemistry 
program^ These points are at every distance R from 
3.25 A to 5.5 A in steps of 0.25 A and for every angle 9 
from 0° to 180° in steps of 15°. This gives a total of 130 
points. We used the RCCSD(T) method (restricted cou- 
pled cluster with single, double and non-iterative triple 
excitations) 31 ! 32 with the aug-cc-pV5Z basis set. 33 i 34 i 35 
The counterpoise procedure of Boys and Bernard^ was 
used to correct for basis set superposition error (BSSE). 
This is the same level of theory and basis set as were 
used in recent calculations of the PES for the A-st&te of 
Ar-SH^ 

Two potential surfaces were obtained from the ab ini- 
tio calculations. These correspond to two adiabatic elec- 
tronic states: one symmetric (A 1 ) and one antisymmetric 
(A") with respect to reflection in the plane of the nuclei. 
The two states are degenerate at linear geometries but 
nondegenerate at nonlinear geometries: the splitting is 
an example of the Renner- Teller effect. The interaction 
energies for each state were interpolated using a 2D spline 
function, and contour plots of the resulting surfaces are 
shown in Fig. ^ The A' and A" surfaces result from 
the electronic Hamiltonian without spin-orbit coupling: 
a discussion of surfaces including spin-orbit coupling is 
presented in section IV. 

The adiabatic surfaces (adiabats) are qualitatively sim- 
ilar to those reported for He-SH and Ne-SH complexes^ 
The latter were calculated at the RCCSD(T) level, using 
the smaller aug-cc-pVTZ basis set, but with an additional 
set of bond functions, and counterpoise correction. A 
comparison of the positions and energies of the minima 
on the Ar-SH surfaces presented here with those for X- 
state He-SH and Ne-SH is given in Table |T| For all the 
A" surfaces there is a global minimum in the linear X- 
S-H configuration (where X is He, Ne or Ar) and a local 
minimum in the linear X-H-S configuration. The X-H-S 
configurations are saddle points on the A' surfaces, which 
have shallow local minima at 6 = 180° and global min- 
ima at nonlinear configurations. The A' global minima 
are deeper than those on the A" surfaces, because in the 
A' state the SH tt hole is directed towards the Ar atom, 
resulting in reduced repulsion. The global minima for the 
A' state occur at angles 9 that increase with the atomic 
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FIG. 1: Ab initio potential energy surface contour plots for 
Ar-SH (X) in the A' state (upper plot) and the A" state 
(lower plot). Solid contour lines are shown at 10 cm -1 inter- 
vals, ranging from cm~' to —150 cm -1 inclusive for the A' 
surface and cm -1 to —120 cm -1 inclusive for the A" sur- 
face. Dashed contour lines are shown at 100 cm -1 intervals 
from +100 cm" 1 to +500 cm -1 inclusive for both surfaces. 
The linear Ar-H-S conformation corresponds to 6 — 0°. 



number of the rare gas atom. Also, as expected, the min- 
ima are deeper for clusters containing heavier (and more 
polarizable) rare gas atoms. 

In order to perform dynamical calculations on Ar-SH, 
we need to evaluate the matrix elements of the poten- 
tial between electronic states labelled with an angular 
momentum quantum number A. For this purpose it is 
convenient to re-express the potential energy surfaces as 
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TABLE I: Positions and well depths of potential minima on 
the A' and A" adiabatic surfaces for the X-state of SH-rare 
gas clusters. The results for Ne-SH and He-SH clusters are 
from ref. H3- 



Cluster 




State 


R/A 


e/° 


depth/cm 1 


Ar-SH 




A' 


3.678 


66.6 


157.69 


Ne-SH 




A' 


3.611 


57.2 


57.05 


He-SH 




A' 


3.639 


54.4 


25.97 


Ar-SH 


A 


and A" 


3.801 


180 


128.54 


Ne-SH 


A 


and A" 


3.593 


180 


54.27 


He-SH 


A 


and A" 


3.593 


180 


25.27 


Ar-SH 




A" 


4.274 





125.22 


Ne-SH 




A" 


4.101 





45.75 


He-SH 




A" 


4.126 





21.16 



the sum (Vb) and difference (V2) potentials 

V (R, 9) = ^ [V A > (R, 0) + Va" (R,e)], 

V 2 (R,e) = ~ [V A > (R, 6) - v A » (R,e)}. 

Contour plots of these surfaces are shown in Fig. They 
are quite similar to those recently presented by Sumiyoshi 
et alM The latter were fitted to reproduce experimental 
results, with starting values for the potential parameters 
obtained from ab initio calculations (RCCSD(T)/aug-cc- 
pVTZ). The form of our sum potential is also qualita- 
tively similar to those recently presented for Ne-SH and 
Kr-SH by Suma et alS^ 

III. BOUND-STATE CALCULATIONS 
A. Coupled channel calculations 
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FIG. 2: Contour plots for the Ar-SH (X) difference (upper 
plot) and sum (lower plot) potential energy surfaces. Solid 
contour lines are shown at 10 cm - intervals, ranging from 
cm -1 to —300 cm -1 inclusive for the difference surface and 
from cm -1 to —120 cm -1 inclusive for the sum surface. 
Dashed contour lines are shown at 100 cm -1 intervals from 
+100 cm -1 to +500 cm -1 for the sum surface only. The linear 
Ar-H-S conformation corresponds to a — 0°. 



The bound states of a complex such as Ar-SH (A 2 n) 
involve both potential energy surfaces. In the present 
work we use a coupled-channel approach to calculate the 
bound states. In a body- fixed axis system the Hamilto- 
nian operator is 

where H mon is the monomer Hamiltonian and V is the 
intermolecular potential. In a full treatment including 
overall rotation, the total wavefunction of the complex 
may be expanded^ 

*l M =R-iJ2*™MX J Pn-,rM(R), (2) 

jPXa 

where the channel basis functions are 



/2J+ i\ 4 
x {-^rj D£ P (a,p,0), (3) 

The monomer basis functions are labelled by Hund's case 
(a) quantum numbers A and er, the projections of the elec- 
tronic orbital and spin angular momentum along the SH 
axis, and u = A+o\ 52 A, a and lo are all signed quantities. 
The D functions are Wigner rotation matrices. 39 The 
first D function describes the rotation of the monomer 
with respect to body-fixed axes, with angular momen- 
tum quantum number j (including electronic orbital and 
spin angular momentum) and projection P along the in- 
termolecular vector R. The second D function describes 
the rotation of the complex as a whole, with total angular 
momentum J and projections M and P onto space-fixed 
and body-fixed axes respectively. The angles (f3, a) de- 
scribe the orientation of the R vector in space. 

The monomer Hamiltonian used here for SH(AT 2 n) is^i 

H moa = bQ -l-s) 2 + H so , (4) 
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where the rotational constant b is 9.465 cm -1 iSi For sim- 
plicity the spin-orbit Hamiltonian H so is taken to be in- 
dependent of R and 9 and equal to aXa, with a = —378.5 

It is convenient to expand the sum and difference po- 
tentials in terms of renormalizcd spherical harmonics 

V {R,6) = J2Vio(R)Cio{d,0); (5) 

i 

V 2 (R,9) = J2V l2 (R)C l2 (6,0). (6) 

i 

The potential matrix elements between the angular basis 
functions may then be written 

(JM;jP; Xa\V\JM;j'P';X'a') 

= 5p P ,5 aa , Vi,\x-y\(R) 3j,a-v jV; P), (7) 
i 

where the potential coupling coefficients are 

gi,\-\< (juj;j'v';P) 

= (-l) p - u [(2j + l)(2/ + l)]» 

The potential matrix elements are independent of J and 
diagonal in P. Nevertheless, in a full treatment the wave- 
functions of the complex are linear combinations of func- 
tions with different values of P, because the operator 
(J — j) 2 in Eq. ^has matrix elements off-diagonal in P 
(AP = ±1). However, the full wavefunctions are eigen- 
functions of the parity operator: symmetrized basis func- 
tions may be constructed by taking even and odd linear 
combinations of 3>/pf\<7 aim ^j-p -x-a- 

In the present work, the coupled equations are solved 
using the BOUND program of Hutson4i The wavefunc- 
tion log-derivative matrix is propagated outwards from a 
boundary point at short range (i? m in) and inwards from 
a boundary point at long range (i? max ) to a matching 
point (i?mid) in the classically allowed region. If E is an 
eigenvalue of the Hamiltonian, the determinant of the dif- 
ference between the two log-derivative matrices at i? m ;d 
is zero42i^ The BOUND program locates eigenvalues 
by searching for zeroes of the lowest eigenvalue of the 
matching determinant^ using bisection followed by the 
secant method. In the present work we use i? m i n = 3.0 A, 
i? m ax = 9.5 A, i?mid = 4.2 A and a log-derivative sector 
size of 0.02 A. The basis set includes all SH functions up 
to jmax = 15/2 in both spin-orbit manifolds. 

The energies obtained from full close-coupling calcula- 
tions for the lowest few J = 3/2 levels of Ar-SH (actu- 
ally carried out in the equivalent space-fixed basis set^) 
are shown in Table |H] These levels all correlate with SH 
2 n 3 / 2 , j = 3/2 and are labelled with the projection quan- 
tum number P and Van der Waals stretching quantum 
number n. We use the convention that levels in which 



TABLE II: Bound-state energies for J = 3/2 levels of Ar-SH 
from full close-coupling calculations (average Ecc and parity 
splitting A£Vjc), helicity decoupling calculations (Ehd) and 
single-surface calculations on the lower adiabatic surface in- 
cluding spin-orbit coupling (Ea<i)- All energies are relative to 
the dissociation energy to form SH (X 2 H 3/ / 2 , j = 3/2). All 
energies are given as wavenumbers in cm -1 . 



p 


n 


Ecc 


A£ C c 




-Ehd 


E A d 


+3/2 





-102.745 


+3.5 x 10" 


-5 


-102.652 


-102.725 


+1/2 





-97.766 


+0.144 




-97.593 


-97.667 


-3/2 





-94.940 


-1.1 x 10" 


-3 


-94.894 


-95.035 


-1/2 





-92.116 


-0.138 




-92.222 


-92.293 


+3/2 


1 


-77.292 


+2.6 x 10" 


-5 


-77.111 


-77.258 


+ 1/2 


1 


-72.276 


+0.134 




-72.100 


-72.148 


-3/2 


1 


-69.356 


-1.3 x 10" 


-3 


-69.265 


-69.382 


-1/2 


1 


-67.207 


-0.124 




-67.291 


-67.315 



P and lo for the dominant basis functions have the same 
sign are labelled with positive P and those where they 
have different sign are labelled with negative Pi In or- 
der of increasing energy, the lowest four levels for Ar-SH 
have P = +3/2, +1/2, -3/2, -1/2, in contrast to Ar- 
OH where the order is +3/2, +1/2, -1/2, -3/2iii The 
difference is due to the anisotropy of the sum potential 
V (R,9): the ratio U20/U10 is larger for Ar-SH. 

The close-coupling results may be compared with 
the microwave experiments of Sumiyoshi et al. ^ who 
obtained a rotational constant B = 1569.66 MHz 
(0.05236 cm -1 ) and parity doubling constant qj = 
0.32873 MHz (1.10 x 10~ 5 cm" 1 ) for the ground state 
(P = +3/2). These correspond to a J = 3/2 — 5/2 sep- 
aration of 0.262 cm -1 and a J = 3/2 parity splitting of 
6.6 x 10~ 5 cm" 1 , which compare with calculated values of 
0.263 cm" 1 and 3.5 x 10~ 5 cm" 1 respectively. The very 
good agreement for the rotational spacing suggests that 
the equilibrium distance of our ab initio potential is quite 
accurate. The difference of almost a factor of 2 in the par- 
ity splitting is less satisfactory, but Dubernet et al~ have 
shown that such terms involve complicated combinations 
of high-order terms involving the difference potential, 
spin uncoupling and Coriolis perturbations. Small dif- 
ferences between the energies of excited states can have 
a large effect on the parity splitting. Sumiyoshi et alm^ 
have very recently measured microwave-millimetre- wave 

double-resonance spectra of the P = +1/2 < h3/2 band 

of Ar-SH. The centre of gravity of the parity components 
of the J — 3/2 <— 1/2 line is 81.8 GHz (2.73 cm- 1 ). The 
corresponding calculated quantity from our potential is 
4.805 cm" 1 . In addition, the measured parity splitting 
for the J = 3/2, P = +1/2 level is about 5300 MHz 
(0.177 cm" 1 ), which compares with 0.144 cm -1 from our 
calculations. An interesting possibility for future work 
would be to adjust the ab initio potential to improve the 
fit to the spectroscopic parameters using the morphing 
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procedure of Meuwly and Hutson4& 

B. Wavefunctions 

The full wavefunctions (Eq. |2J) contain contributions 
from all possible values of P and uj and are not separable 
between the body-fixed angles (6, 4>) and the space-fixed 
angles (/3, a). This makes them hard to visualize. In 



addition, since the mixings depend on the total angu- 
lar momentum J, they are not convenient for calculating 
band intensities. We therefore introduce two approxi- 
mations to simplify the description of the wavefunctions 
for this purpose. First, we introduce the helicity decou- 
pling approximation, where matrix elements of (J — j) 2 
off-diagonal in P are neglected. Secondly, we neglect ma- 
trix elements of H mon off-diagonal in a (spin-uncoupling 
terms). The coupled equations then simplify to 



ti 1 d 2 
2/i dR 2 



771111011 



2fiR 2 



(J(J+1)+3(3 + 1)-2uj 2 )-E1 



= -J2 (JM;j'PX'a\ V \JM;jPXa) x£„ ;i ,. AV (i?). 
j>\> 



(9) 



Since all matrix elements off-diagonal in P and a have 
been neglected, states with quantum numbers (P, uj) and 
(— P, —uj) arc uncoupled and it is not necessary to take 
combinations of definite total parity. However, states 
with (P, uj) and (-P, uj) or (P, —uj) have different poten- 
tial energies and are nondegenerate. 

The energy levels obtained from helicity decoupling 
calculations for Ar-SH are included in Table [H] The 
approximation is accurate to about 0.2 cm -1 for n = 
and 1, but is less reliable for higher states. In particu- 
lar, the region between —60 and —40 cm -1 contains both 
j = 3/2, n = 2 and j = 5/2, n — levels. In the pres- 
ence of the resulting near-degeneracies, the terms that 
are omitted in the approximate Hamiltonian can cause 
quite significant level shifts. 

In the helicity decoupling approximation, P is a good 
quantum number. However, uj is not because V% mixes 
levels with AA = ±2 (but Ac = 0) and thus mixes uj — 
+3/2 with -1/2 and uj = +1/2 with -3/2. However, 
in the absence of terms off-diagonal in a the two sets are 
not mixed with one another. Each wavefunction thus has 
only two components corresponding to different values of 
uj. The wavefunctions may be written 



(10) 



where the basis functions now exclude the 0-dependence, 
/2J+1^ /2 



- ip a <p X 



V 8tt 2 



DjSpifx, /3,<f>), (11) 



and the 2-dimensional functions that characterize the 
components of the wavefunction for each uj are 

XPn;JR,0) = £ (j + 5) 1/2 4j0)XP»y;Aa(#), (12) 



where cPj U (9) is a reduced rotation matrix. 
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We have adapted the BOUND program 41 to calcu- 
late wavefunctions for this case by back-substituting into 
the log-derivative propagation equations, as described 
for the closed-shell (single surface) case by Thornley 
and Hutson^ 7 Examples of the resulting wavefunctions 
Xp„. w (fi, 9) are shown in Fig. EI It m &y be seen that the 
components for different values of uj have quite different 
radial and angular distributions. 

The Ar-SH wavefunctions are qualitatively similar to 
those for Ne-SH obtained by Cybulski et al& Since the 
potential anisotropy for Ar-SH is only a few tens of cm -1 
in the well region, there is only weak mixing of SH ro- 
tational functions with different values of j. For this 
reason the wavefunctions are dominated by the functions 
dp^(6), as described by Dubernet and Hutsoni for the 
case of Ar-OH. The d functions are shown, for example, 
111 Fig. 7 of ref.Q and the angular parts of the wavefunc- 
tions of Fig. |21 follow them quite closely. 



C. Adiabatic approximations 

In a basis set of Hund's case (a) functions with signed 
values of A = ±1 and a = ±1/2, we can define new adia- 
batic surfaces (adiabats) including spin-orbit coupling as 
eigenvalues of the Hamiltonian matrix at each value of R 
and 6, 



( 



\ 



0^2' 


V-2 









Vt 



- 5 ( 




V 



Vi 


' 





v 2 



V + ^aJ 



(13) 



where again the spin-orbit Hamiltonian is taken to be 
simply aXa. This clearly factorizes into two equivalent 
2x2 matrices, one containing uj = +3/2 and —1/2 and 
the other containing uj = +1/2 and —3/2. The resulting 
adiabats may be designated V+(R, 9) and V-(R, 9) with 
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30 60 90 120 150 180 30 60 90 120 150 180 
/degrees G /degrees 

FIG. 3: [Colour Figure] Contour plots of the wavefunction components, superimposed on the average potential. 



corresponding electronic functions (R, 9) and ip- (R, I 
given by 



cosa ad (i?,6 
- sin a ad (R, 



sin a ad {R,t 
cos a ad (i?,/ 



^±3/2 n 4) 



where ip u = ip\tpo-. The adiabats for Ar-SH are shown 
in Fig.0] Since for Ar-SH V2(R, 9) is small compared to 
a in the well region, the lower adiabat is always predom- 
inantly lu = ±3/2 in character and the upper adiabat 
is always predominantly to — =Fl/2 in character. The 
corresponding mixing angle a a( j is at 9 — and 180° 
(where Vz(R, 9) — 0) and less than 20° at other angles for 
i? > 3.5 A. A contour plot of the mixing angle is shown 
in Fig. 

A further consequence of the large spin-orbit coupling 
constant is that both adiabats resemble the sum potential 
Vq(R, 9) much more than the A 1 and A" potentials. The 
spin-orbit coupling has in effect quenched the splitting 
between the A' and A" states. This explains why there 
is no tendency for the wavefunctions shown in Fig. [31 to 
"fall into" the non-linear minimum of the A' state. 

The existence of adiabats including spin-orbit coupling 
suggests a Born-Oppenheimer separation in which the 
total wavefunctions are written approximately as 



R-%{R,e) X iPn(R,6), 



(15) 



where ipi is one of the functions of Eq. PH] and Xin (R, 9) 
is a solution of an effective Schrodinger equation of the 
form 



ti z d 2 

2fi dR 2 



Vi(R, 9) — Eip n 



2^R 2 



Xi Pn(R,0)=0. 



(16) 



However, the appropriate angular operator H rot to use in 
such a calculation is hard to define. The reduced rotation 
matrices d 3 Puj (9) that describe the free SH molecule are 
eigenf unctions of H lot — b(j 2 — 2uj 2 ), where 



r 



d 



sin 9 89 



89 



P 2 



2Plo cos 9 



(17) 

This contains singularities at 9 — and/or 180° that de- 
pend on the values of P and u>. However, there is no 
single value of oj that is appropriate at all configurations. 
The simplest approach is to replace to in Eq. El with the 
value that is appropriate at 6 = and 180° and solve 
Eq. in a basis set of d functions for each value of P. 
This is equivalent to solving the coupled equations using 
a basis set containing only functions with a single value 
of u). The results obtained with this approximation are 
included in Table [H] it may be seen that it gives energies 
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FIG. 4: Contour plots for Ar-SH (X) adiabats including spin- 
orbit coupling. Contour lines are shown at 10 cm - intervals, 
ranging from 70 cm -1 to +a/2 for the upper surface and from 
—310 cm" 1 to —a/2 for the lower surface. 



that are generally slightly too low (compared to the helic- 
ity decoupling results), by 0.05 to 0.15 cm" 1 . A slightly 
better but significantly more complicated approximation 
would be to replace to with (u>) and oj 2 with (lu 2 ) in Eq. 
1171 to give an improved effective potential. 

One approach that is clearly not appropriate is to 
carry out a bound-state calculation on a single adiabat 
V± (R, 9) assuming that the SH molecule behaves as a 
closed-shell rigid rotor. Such a calculation would give 
substantially incorrect energies and wavefunctions. 

It is in fact true that no wavefunction of the form \ 15\) 
can have the correct behavior near both linear geometries. 
To see this, consider an alternative definition of the mix- 
ing angle that can be obtained from a single wavefunction 



8 



4.0 



< 

DC 



3.5 




60 90 120 
/ degrees 



180 



FIG. 5: Contour plot of the adiabatic mixing angle a ac j • This 
angle is derived from the adiabats and is defined in Eq. If 41 
The contour lines are spaced at 2° intervals. 



in the helicity decoupling approximation, 
tana J Pn (R,9) = — 



X 



Pn;±3/2 



(R, 



(18) 



This quantity is plotted for n — and all four P values 
corresponding to j — 3/2 in Fig. El The mixing angles 
for P = +3/2 and P = +1/2 bear some similarity to 
a a d(i?, #) (Fig. |SJ at small 9, but tend to 90° instead 
of zero at 9 = 180°. Conversely, the mixing angles for 
P = -3/2 and P = -1/2 tend to 90° at 9 = 0. This 
is easy to explain in terms of the reduced rotation ma- 
trices that appear in Eq. 1121 For example, the functions 
^±3/2 ±3/2 W au behave as cos 3 (9/2) as 9 — > 180°, while 
the functions d 3 ^ /2 +a/2 (9) behave as cos(0/2). This cor- 

90° 



ooasfl^ 180° so a 



J 

3/2, n 



=Fl/2,±3/2 

responds to tan a^, 2 n 
in that limit. The point here is that the component of 
the P = +3/2 wavefunction on the u> — —1/2 surface 
goes to zero more slowly than that on the u = +3/2 
surface as 9 —> 180°. Fig. shows that the coupled- 
channel wavefunctions (|10fl for P — +3/2 and +1/2 are 
predominantly in the u> = —1/2 state near 9 = 180°, 
which corresponds to the upper adiabat rather than the 
lower one. The P = —3/2 and —1/2 wavefunctions show 
similar behavior around 9 = 0. This is not the behavior 
implied by Eq. 



IV. ELECTRONIC SPECTRUM CALCULATION 

A. Transition wavenumbers 

In order to calculate the line positions in the vibra- 
tionally resolved A 2 E + <— X 2 II electronic spectrum, 
we require the bound-state energies for the excited elec- 
tronic state, as well as those for the ground state. For 
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FIG. 6: [Colour Figure] Comparison of the adiabatic mixing 
angle a a d (thick black line) with angles obtained from wave- 
functions correlating with j — 3/2, u — +3/2 for P = +3/2 
(red), P = +1/2 (blue), P = -1/2 (green) and P = -3/2 
(black). The mixing angles are shown as cuts through the 
corresponding surfaces taken at R — 3.9 A. 



the A 2 S + state we make use of the recent PES presented 
by Hirst et alm^ This surface has a global minimum of 
-742.5 cm -1 at the linear Ar-H-S conformation (9 = 0°), 
and a secondary minimum of —673.7 cm -1 for linear Ar- 
S-H [9 — 180°). The two minima are separated by a 
barrier more than 600 cm -1 high and the lowest-energy 
vibrational levels are localised within one or other of the 
two wells. 

The bound states of this PES have been analysed 
previously^ and only a brief discussion is given here. 
Bound-state energies were calculated as eigenvalues of 
the spin-free triatomic Hamiltonian in Jacobi coordi- 
nates. Discrete variable representations (DVRs) were 
employed for both the intermolecular distance R and the 
angle 9. For R, 120 sinc-DVR functions^ were used, with 
DVR points ranging from 2.5 A to 8.5 A. For 9, a 64-point 
DVR based on Legendre polynomials was used. With this 
basis set, the bound levels of interest were converged to at 
least seven significant figures. The resulting levels are la- 
beled by quantum numbers (i>sh, b K ,n), where i>sh an d fi 
are quantum numbers for the SH stretch and the atom- 
diatom stretch respectively. K is the projection of the 
total angular momentum of the diatom, neglecting spin, 
onto the body-fixed z-axis, and b is the number of nodes 
in the intermolecular angle 9. The resulting energies for 
levels with total angular momentum N — (neglecting 
spin) are in precise agreement with previous results^ To 
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facilitate the calculation of band intensities, for N = 1 
the helicity decoupled approximation was employed, in 
which the Coriolis terms coupling different i^-levels are 
ignored. The helicity-decoupled energies are within 0.5 
cm -1 of the full close-coupled results^ 

For the purpose of calculating transition frequencies, 
the asymptotic separation of the potentials is taken to be 
the experimental excitation energy from the v = 0, j = 
3/2 level of the 2 n 3 / 2 state to the lowest v = 0, j = 1/2 
level of the 2 £+ state of isolated SH, which is 30832.68 
cm -1 4& All transitions of the complex were assumed to 
originate from the P = +3/2 level of Ar-SH ( 2 II). The 
lowest-energy transition frequency for the complex is cal- 
culated to be 30488.5 cm -1 , which is 31.5 cm -1 greater 
than the experimental value of 30457 cm^ 1 .^ This agree- 
ment is reasonable, considering the level of theory used 
in the calculation of the potentials. 



B. Transition dipole moments 

Calculations of spectroscopic intensities require tran- 
sition dipole moments fj.f ot , where 

rf ot = (i|/i e l|f) . (19) 

The integrals involve the initial (i) and final (f) wave- 
functions as determined from bound-state calculations. 
In this work we evaluated transition dipoles over internal 
coordinates (R, 9), neglecting overall rotation. This gives 
transition dipoles that correspond to band intensities be- 
tween intermolecular vibrational states. The electronic 
dipole moment, fi e i, is in general a parametric function 
of the nuclear coordinates. In the body-fixed frame it 
may be expanded in terms of reduced rotation matrices, 

^{R,9) =Y^^( R ) d AP,AxW> ( 20 ) 

j 

where AP = P\ — Pf and AA = Ai — Af. In the ex- 
cited electronic state, P = K ± |. In this work it is 
assumed that /i c i consists purely of contributions from 
the SH monomer, so that only j ' = 1 contributes in Eq. 
1201 and the coefficients are independent of R. Since 
we are dealing with a perpendicular transition in SH, 
AA = ±1. The transition dipoles were calculated as one- 
dimensional Gaussian quadratures in 9, then integrated 
over R. 



C. Intensities and lifetime factors 

The signals in a pulsed-laser fluorescence excitation ex- 
periment decay exponentially following each pulse, with 
a lifetime equal to that of the excited state being probed. 
Intensities are typically measured as the area under the 
decay curve, and so the experimental intensities are pro- 
portional to the lifetime of the excited state. Ar-SH is 



somewhat unusual in the large range of lifetimes exhib- 
ited by different vibrational levels in the A state. It 
is known that the presence of the Ar atom blocks the 
electronic predissociation of the SH radical, leading to a 
greatly increased lifetime of up to 600 ns for low-lying 
bound levels, compared to ~ 1 ns for the uncomplexed 
species i^S^ However, the actual lifetime depends on the 
degree of vibrational excitation, and lifetimes specific to 
particular levels have been calculated by McCoy^S 

The top panel of Fig. [7] shows a spectrum calculated 
directly from the squares of transition dipoles, while the 
center panel shows a spectrum in which the intensities 
have been multiplied by McCoy's lifetime values. Clearly 
this is possible only for levels for which lifetime data exist, 
and transitions to other levels are omitted in the center 
panel (i.e., it is assumed that their lifetimes are small). 
Also shown is an experimental spectrum from Ref. [2(J. 
The spacings between the peaks in the calculated Ar-SH 
spectrum are consistently ~ 5% smaller than in experi- 
ment. It is clear that the intensity distributions are sig- 
nificantly different in the two calculated spectra, and that 
the one that includes lifetime factors gives considerably 
better agreement with experiment. The agreement in 
intensities is quite good, especially considering that the 
experimental spectrum was most likely not normalised 
for dye laser power^i From our results it seems likely 
that the small peak at ~ 30810 cm" 1 in the experimental 
spectrum can be assigned to the transition to (0, 0°, 6). 

Even without the lifetime weighting, transitions to lev- 
els localized in the 9 — 180° well of the A state, which 
are labelled with + symbols in Fig. [7| are weak in the 
simulated spectrum. This arises because of poor overlap 
with the P — +3/2 ground-state wavefunction (which is 
concentrated around 9 — 0). In addition, it is likely that 
such levels have short lifetimes close to that of uncom- 
plexed SH22i and so will have even lower intensities in the 
fluorescence excitation spectrum. These levels have not 
been observed experimentally to our knowledge. 



V. SUMMARY 

We have obtained new ab initio potential energy sur- 
faces for the Ar-SH complex in its ground 2 n electronic 
state and used them to calculate bound-state energies 
and wavefunctions using coupled-channel methods. We 
have also described a new adiabatic approximation that 
includes spin-orbit coupling and can be used to calculate 
bound states on a single potential energy surface. How- 
ever, the adiabatic wavefunctions fail to reproduce some 
features of the true wavefunctions. We have used our re- 
sults to simulate the vibrationally resolved laser-induced 
fluorescence excitation spectrum of Ar-SH, with intensi- 
ties modelled using calculated transition dipole moments 
and calculated lifetimes. The inclusion of the lifetime 
factor is important to obtain satisfactory agreement with 
the experimental intensities. 
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FIG. 7: Calculated vibrationally resolved fluorescence- 
excitation spectrum of Ar-SH for the A 2 E + - X 2 II electronic 
transition. The lifetime weighting of the intensities is absent 
for the the top panel, and present for the middle panel. Lines 
labeled with + indicate transitions to the 9 — 180° well of the 
A state. For comparison, the experimental spectrum is shown 
in the bottom panel, taken from Ref. l2pl Note that the ex- 
perimental spectrum contains contributions from Ar2-SH and 
larger clusters as well as Ar-SH. 
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